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Inelastic tunneling effects on noise properties of molecular junctions 

Michael Galperin,^ Abraham Nitzan,^ and Mark A. Ratner^ 

^Department of Chemistry and Nanotechnology Center, Northwestern University, Evanston IL 60208 
'^School of Chemistry, The Sackler Faculty of Science, Tel Aviv University, Tel Aviv 69978, Israel 

(Dated: February 6, 2008) 

The effect of electron-phonon coupling on the current noise in a molecular junction is investigated 
within a simple model. The model comprises a 1-level bridge representing a molecular level that 
connects between two free electron reservoirs and is coupled to a vibrational degree of freedom 
representing a molecular vibrational mode. The latter in turn is coupled to a phonon bath that 
represents the thermal environment. We focus on the zero frequency noise spectrum and study the 
changes in its behavior under weak and strong electron-phonon interactions. In the weak coupling 
regime we find that the noise amplitude can increase or decrease as a result of opening of an inelastic 
channel, depending on distance from resonance and on junction asymmetry. In particular the relative 
Fano factor decreases with increasing off resonance distance and junction asymmetry. For resonant 
inelastic tunneling with strong electron-phonon coupling the differential noise spectrum can show 
phonon sidebands in addition to a central feature. Such sidebands can be observed when displaying 
the noise against the source-drain voltage, but not in noise vs. gate voltage plots obtained at low 
source-drain bias. A striking crossover of the central feature from double to single peak is found 
for increasing asymmetry in the molecule-leads coupling or increasing electron-phonon interaction. 
These variations provide a potential diagnostic tool. A possible use of noise data from scanning 
tunneling microscopy experiments for estimating the magnitude of the electron-phonon interaction 
on the bridge is proposed. 

PACS numbers: 72.70.-|-m 72.10.Di 73.63.-b 85.65. -|-h 
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I. INTRODUCTION 



Development of molecular electronics as a complement 
to traditional semiconductor-based electronics provides 
challenges both experimentally and theoreticallyiSiSi^ 
Early studies of molecular junctions were restricted to 
measurement of current- voltage (or conductance- voltage) 
characteristics of such junctions and their dependence on 
junction parameters such as wire length, molecular struc- 
ture, molecule-lead coupling, and tcmpcratureALSiSiiSiii 
More recently, charging (Coulomb blockade) and inelas- 
tic and mechanical effects appear at the forefront of 
researchiiSiiiiiiiSiiSiii In particular, inelastic electron 
tunneling spectroscopy (lETS) has become an important 
diagnostic and characterization tool, indicating the pres- 
ence and structure of a molecule in the junctioniSiii and 
providing information on the junction structured 

In addition to the I/V behavior, noise char- 
acteristics can provide information about molecular 
junctions^ Noise measurements in mesoscopic tun- 
nel junctions have been under study for a long fine, 
e.g. in semiconductor double barrier resonant tunnel- 
ing structures (DBRTS))^":^^'^^'^'^ quantum point con- 
tacts (QPC)3i:^5:^Mi Coulomb blockaded Josephson 
junctions 



28.29 



and advances in understanding their di- 
agnostic capabilities are being made (see, e.g. b(J). Re- 
sults of measurements in molecular tunnel junctions are 
expected to appear in the near future, though these 
measurements are difficult due to the background 1// 
noise caused by charge fluctuations in the environ- 
ment. Many theoretical studies of noise in nano junc- 
tion transport are available in the literature (see Ref. ili^ 



and references therein). In particular noise in DBRTS 
was studied within the nonequilibrium Green function 
(NEGF) formalism in the ballistic'^^-'^^ and the Coulomb 
blockadaSiM regimes. Noise in QPC was studied within 
NEGF^^ and kinetic equation approachespS, while noise 
in Coulomb blockaded Josephson junctions was inves- 
tigated within a phase correlation theory approachi^ 
NEGF was also applied to study shot noise in chain 
models'^*'^^ and disordered junctions.""^ Shot noise for co- 
herent electron transport in molecular devices was inves- 
tigated within a scattering theory approach in a num- 
ber of works 4i*^2ii^ Inelastic effects in the noise spec- 
trum were studied recently, first in connection with nano- 
electromechanical systems (NEMS)^^*^^i^S4LiS4^ and ac- 
driven junctions^Si^ Substantial work on this issue has 
been done within a scattering theory approach pSiSi^iS 
Finally predictions on phonon assisted resonant shot 
nois e sp ectra of molecular junctions were obtained in 
Ref. Hif within NEGF. 

In this work we investigate inelastic tunneling and its 
influence on the zero frequency noise in a simple one- 
level molecular junction model. The one-level junction 
model focuses on one molecular orbital (say the lowest 
unoccupied molecular orbital, LUMO) that supports res- 
onance transmission beyond a certain voltage bias thresh- 
old. The molecular orbital is coupled to a local oscillator 
representing the molecular nuclear subsystem. The oscil- 
lator in turn is coupled to a thermal harmonic bath. We 
have recently used this model within the NEGF method- 
ology to study inelastic effects on electronic transport 
in molecular junctions in the weak^^ and in the inter- 
mediate to strong— electron phonon coupling regimes. 
Here we utilize the same approaches to investigate in- 



elastic effects in the noise spectrum. The physics of the 
problem is determined by the relative magnitudes of the 
relevant energy parameters: A_E - the spacing between 
the leads' Fermi energies and the energy of the molecu- 
lar orbital, F - the broadening of the molecular level due 
to electron transfer interaction with the leads, M - the 
electron-phonon coupling and lo - the phonon frequency. 
We consider both resonant and far off-resonant inelastic 
tunneling cases. In the off-resonant limit, /S.E ^ F it 
is usually the case that also /^E ^ M and the process 
can be described in the weak electron-phonon coupling 
limit. For resonance transmission, AE < F, the electron- 
phonon interaction may still be weak, M ^ F^fiS These 
weak coupling situations can be treated within standard 
diagrammatic perturbation theory on the Keldysh con- 
tour, see e.g. 1571 More interesting physics is observed in 
the case of strong electron-phonon interaction, M ^ F, 
which we treat within the self-consistent NEGF scheme 
introduced in Ref. t58t In both limits we apply the cor- 
responding approach to study effects of electron-phonon 
coupling on the zero- frequency noise and the Fano factor. 

Our results are discussed in comparison with recent 
work on the effect of electron-phonon coupling on the 
noise properties of molecular junctions. In the weak cou- 
pling limit our result agree qualitatively with the scat- 
tering theory based treatment of Chen and Di VentraS^ 
in the case of near resonance transmission in symmetric 
junctions. However, in contrast to Ref. |5a that finds that 
the leading contribution to this effect is of order M^ we 
find that this effect is of order AP. Also, different qualita- 
tive behavior is found in off-resonance situations and for 
resonance transmission through asymmetric junctions. 

In the strong electron-phonon coupling limit we com- 
pare our result to that of Zhu and Balatsky^^ who treat 
this problem using the (essentially scattering theory) pro- 
cedure of Lundin and McKenzie>^ Again we find some 
new qualitative effects that can not be obtained from the 
simpler approach. 

Our analysis suggests the usefulness of the "noise spec- 
troscopy" , either in a source-drain response or in a gate 
voltage response, to measure the values of the electron- 
vibration coupling on the bridging molecule and to char- 
acterize the asymmetry in the coupling between the 
molecule and the source and drain electrodes. 

Section|n]describes formal derivation of the expression 
for the noise through a molecular junction. In Section Hill 
we introduce the model. Section IIVI briefiy reviews the 
noise properties of the elastic current. SectionlVlsketches 
a procedure and presents numerical results for the case of 
weak electron-phonon coupling and Section IVTl does the 
same for the strong coupling case. Section lVlIl concludcs. 



II. GENERAL NOISE EXPRESSION 

A general expression for the frequency dependent noise 
in a 2-terminal junction can be obtained by applying the 
method of Meir and Wingreen. In its most general formu- 



lation this method assumes a single electron exchange be- 
tween bridge and leads and invokes the non-crossing ap- 
proximation (NCA)Si for the bridge-leads coupling. The 
system is described by the Hamiltonian 



H — Hm + Hl + Hji 



E E 



t^/crnC^C^TT, + xl.C. 



(1) 

where Hm, Hl and Hn denote the Hamiltonian of the 
independent molecule and the left and right leads, re- 
spectively and where Cm (cj„) and Ck (cj.) denote sin- 
gle electron annihilation (creation) operators in molecule 
and leads. The current operators for the molecule and 
each lead are defined by 

iK(i) = -^E E {ykmcl(t)c.^{t)-Vmkc\^{t)Ck{t) 
k&K m£M 

(2) 
and the average steady state current for a given potential 
bias is obtained in the form.®^-^'^ 



'A' 



'K 



dE 

2^ 



ik{E) 



(3) 



2e^ 



^K{E) = -Tr[S<(£;)G>(£;)-S>(i?)G<(i?)](4) 

where S^'^ are lesser/greater projections of the self- 
energy due to coupling to the (metallic) contact K [K = 
L,R) and G^'^ are lesser/greater Green functions, all 
these are defined in a molecular subspace of the prob- 
lem. Obviously II = —Ir in steady state. 

The noise spectrum is defined as the Fourier transform 
of the current-current correlation function. The relevant 
current fiuctuations have to be considered carefully be- 
cause by their nature such fiuctuations correspond to 
transient deviations from steady states. We follow the 
analysis of Ref. JGJ that leads to the following expression 
for the operator that corresponds to the current in the 
outer circuit 



Ht) = vMt) + vnMt) 



(5) 



where 



Vl = 



Cr 
C 



riR = 



Cl_ 
C 



C = Cl + Cr (6) 



Cl and Cr are junction capacitance parameters that 
take into account the charge accumulation at the corre- 
sponding bridge-lead interface. Note that at steady state 
I = II = —Ir. The different signs assigned to rjL and tjr 
stem from the fact that the current is considered positive 
on either side when carriers enter the system (molecule) . 
In terms of /, Eq. (jHJ, the noise spectrum is defined by 



+ 00 

2 / J-f au\„i'^t 



S{uj) = 2 / dtS{t)e 

J — oo 

Sit) = i({A/(t);A/(0)} 



(7) 
(8) 



where A/ = 1 — 1 and where the averaging is done as 
usual (and as in Eq. ||2Jl) on the non- interacting state of 
the system at the infinite past. The same approacliSSiSi 



J 



that leads to Eq. Q now leads to the following general 
expression for the noise 



Ki,K2 = {L,R} 



/+00 irp 



G<(E- uj) pK, GEkJ> (E) + G>{E~ uj) [E^, GEkJ< (E) 

[J:k, Gf{E- u;) [J:k, G]> (E) - [GEkJ< {E - Lj) [GEkJ> {E) + {lu^ -uj)] 



(9) 



A similar approach was used by Bo and Galperin®^ to obtain an expression analogous to © for a single level bridge 
modeli^ In this paper we restrict our consideration to the case of zero frequency noise which is the relevant observable 
when the measurement time is long relative to the electron transfer time, a common situation in usual experimental 
setup. Then the expression simplifies to 



S{uj = 0) 



4e" 



2 /•+00 



dE 

2^ 



J2 {m.m.Tr [Sk,,k. [G< S>^ + E<^ G>] - 2E<^ ^^JG^\']) ~ i\E) 

Ki,K2 = {L,R} 



(10) 



where i{E) = YjK={l.,r} VKixiE) with ixiE) given by Eq. IJU. 



III. THE SINGLE LEVEL BRIDGE/FREE 
ELECTRON LEADS MODEL 



Eq. (|1U|) can be applied, under the approximations 
specified to very general bridge and leads models. Fur- 
ther progress can be made by specifying to simpler mod- 
els. As usual we assume that the contacts are treated 
as reservoirs of free electrons, each in its own equilib- 
rium. The molecule is described by the simplest coupled 
electron-phonon system as follows. Assuming that the 
energy gap between the orbitals relevant for a molec- 
ular junction (in particular the HOMO-LUMO gap) is 
much greater than level broadening due to coupling to 
the contacts, we consider the electron current through 
each orbital separately. Thus the molecular junction is 
represented by a single level (molecular orbital) coupled 
to two contacts (L and R). The electrons on the bridge 
are coupled to a local normal mode (henceforth referred 
to as primary phonon) represented by a harmonic os- 
cillator. This is coupled to a thermal bath described 
by a set of independent harmonic oscillators ( "secondary 
phonons" ) . The Hamiltonian is thus given by 

H = eoc'c+ ^ ekclck+ ^ fVfcC^c + h.c. 

ke{L,R} ke{L,R} 

+ hujQO^a -I- ^ TlUJpVpbp 4- MaQaC'C + ^ UpQaQp 

p p 

(11) 

where & (c) are creation (destruction) operators of elec- 
trons on the level, c|, (c^) are the corresponding operators 
for electronic states in the contacts, (V (a) are creation 



(destruction) operators for the primary phonon, and &! 
{hp) are corresponding operators for phonon states in 
thermal (phonon) bath. Qa and Qp are shift operators 



^t 



Qp 



5/3- 



fl2) 



Detailed discussion of the model and our self-consistent 
calculation procedures are presented in Refs. |53 andlSSi 
for weak and strong electron-phonon couplings respec- 
tively. Sections Ivl and IVTl below provide brief outlines of 
these procedures. 

In both weak and strong coupling cases the calcula- 
tion is aimed to obtain the electron Green function (or 
rather its Langreth projections^^ on the real time axis). 
Those Green functions are used for the calculation of the 
steady-state current across the junction ^ and the noise 
spectrum (|1(JII . The general current and noise expression 
remain as before, Eqs. (0), (jHl and H10() (excluding the 
trace operations). The lesser and greater self-energies as- 
sociated with the bridge-contact coupling are now given 
by 

Y,<{E) = ijK{E)VK{E) (13) 

S>(i?) = -^[l-fK{E)]VK{E) (14) 

with fx^E) the Fermi distribution in the contact K = 
L, R and 



Tk{E) = 2^Y. \Vk?5{E-ek) 



(15) 



keK 



In the calculation described below we adopt the wide 
band approximation in which T k is assumed constant. 
In this approximation EJf=[E^] = — ir/f/2. 



The observed conduction and noise properties of our 
system depend on parameters that can not be accounted 
for by the present single-electron level of treatment. The 
capacitances Cl and Cr that define the parameters tjl 
and "qR that enter the noise calculation belong to this 
group. We also use below the voltage division parame- 
ter 5 that determines the voltage induced shifts in the 
leads' electrochemical potentials ijll and /x/j relative to 
So according to 



^iL = Ep + SeV = Ef + cVl 

/ifl = Ef - {I - 5) eV = Ef + eVR 



(16) 



This parameter affects the way by which the externally 
imposed bias translates into the relative positioning of 
electronic energies. 

In addition, consider the ratio 



the source-drain voltage V = {^jll — IJ-r)/^ is kept small 
and fixed, fc^T ~ eV ^ F, and the gate voltage is varied 
in order to move the molecular orbital energy in and out 
of resonance. In the second, the gate voltage Vg is fixed, 
e.g. Vg = 0, and V varies so that resonance conditions 
are usually achieved for ksT <C F ^ eV . In the latter 
case we restrict our consideration to the case where the 
onset of resonance transmission is caused by ^l crossing 
the LUMO for a positively biased junction, and disre- 
gard other possible scenarios such as the jir crossing the 
HOMO.^^ In both cases taking ksT ^V provides better 
resolution of resonance features. Also, in the resonant 
tunneling situation it is convenient to consider the first 
derivatives of the current and noise with respect to volt- 
age, dl/dV and dS{uj = 0)/dV respectively, rather than 
/ and S{uj = 0) themselves. 



Tl ^ Tr 

a — — 1 — a = — 
F F 



(17) 



IV. NOISE IN ELASTIC TRANSPORT 



that represents the asymmetry in molecule-lead couplings 
(wide band limit is assumed here and below). This 
is an experimentally controlled parameter that can be 
changed, e.g. with tip-molecule distance in an STM con- 
figuration or by using thin insulating layers to separate 
molecule from leads iS£^ 

The connection (if any) between the parameters a, 6 
and T] = f]L {t]r = T]L — 1) is obviously of great interest. 
An equivalent circuit model for a molecular (or nanodot) 



junctior 



,64.69.70 



usually describes the junction as serially 



connected RC circuits. Assuming that at steady state 
charges ql and qr accumulate on the corresponding ca- 
pacitors the steady state relationships Vl = Ql/Cl — 
IR^ = TLqLRL imply that Tl = {RLCLy\ These re- 
lations imply in turn that 



1-S 1 



V 



1 



V 



(18) 



This leaves two important parameters in the model. 
These however can, at least in principle be determined 
experimentally: The rates F can be determined by op- 
tical experiment*^"'^ and the potential distribution on a 
conducting junction can be probed as demonstrated in 
Ref. 172. In the present paper we use a = F^/F to char- 
acterize the junction asymmetry and several arbitrary 
values oi T] — rjL — Cl/C chosen to represent different 
possible manifestations of noise behavior. 

In addition, in this paper we discuss two experimental 
manifestations of nanojunction transport. In the first. 



Before discussing inelastic effects we briefly review the 
issue of noise in the elastic tunneling situation. In this 
case current, Eq.Q, reduces to the well known Landauer 
expression 



+00 



dE 

2^ 



n{E)[fL{E)-jR{E)] (19) 



while the zero frequency noise expression l(Tn|l reduces to 
a sum of thermal contribution, St = St{uJ = 0), due to 
thermal excitations in the contacts and shot noise term, 
Ss = 5s(cj = 0), due to the discrete nature of the electron 



transporti^. 



S 
St 

Ss 



= S{uj = 0) = St{Lo = 0) + Ss{lu = 0) (20) 

X {fL{E)[l - hiE)] + fR{E)[l - fR{E)]} 
4e2 r+°° dE 



(21) 



2tt 



ToiE) {1 - ro(^)} 



(22) 



with 



X [fLiE)-fRiE)r 



To{E) 



TlTr 



{E-eoy + (T/2r 
Using Eq. ^ this leads to (with /^^^ = fc^T) 



(23) 



J 



dS^ 
dV 

dSs 
dV 



4e3/3 

h 






. sinh(/?(^ - ^iL)/2) 
coah\l3{E~^iL)/2) 



-il-S) 



sinh(/3(^ - iir)/2) 



cosh%(3{E-fiR)/2) 



^°° ^ ToiE)[l niE)] 2[U{E) - fRiE)] 



1 



4cosh^(/3(-B - a*l)/2) 4cosh^(/3(-B - ^IR)/2) 



(24) 
(25) 



Interestingly, these results do not depend on the capaci- 
tance ratios rj. These expressions can be further simpli- 
fied in two cases: (a) When the source-drain voltage V 
is varied we consider, as indicated above, the situation 
where the signal is dominated by the resonance associ- 
ated with the crossing between eg and ^l- In this case 
the terms involving ^jf in Eqs. H24|) and H25|l are disre- 
garded, (b) When the gate voltage Vg is varied at small 
V the two terms on the r.h.s. of Eqs. H24|) and H25|l are 
essentially the same and can be combined together. 

In these elastic transmission cases both experimen- 
tal procedures discussed above give essentially the same 
information. Recall that conductance dl/dV is a 
Lorentzian function of V that peaks at the position of 
the resonance level, with a width (FWHM) related to 
the total escape rate F. In contrast, the differential 
noise graph {dS/dV vs. V) depends on the asymme- 
try of the molecule coupling in the junction a. St itself 
is a Lorentzian function of V, so dSt/dV vs. V has the 
form of a Lorentzian derivative that vanishes at the level 
position and peaks at Eq T r/2 with peak/dip heights 
^ ±^-j^^^Y^ respectively i2^ The contribution from the 
zero temperature shot noise, dSg/dV vs. V , depends 
on the coupling asymmetry a. A double peak structure 
appears when Tq{E)[1 — T[){E)] has three extrema as a 
function of E. The condition for this to occur is 



eV (eV) 
4.0 




i<o 



(26) 



-O.I O.I 

/iL-eo (eV) 



so that for 



l-l/\/2 



< a < 



1 + 1/V2 



it yields a two peak 



structure with peaks at i? = Eq T [r/2]-\/8Q!(l — a) — 1 
and peak heights of 1/4 (in units of 4e^/?i). This behavior 
of ^^(T^) is associated with the transition from double- 
barrier resonance structure at moderate asymmetry to a 
single tunneling barrier for the strongly asymmetric case. 
Ss vanishes when Tg = or Tq = 1 and maximizes when 
To — 1/2. This leads to a two peak structure in sym- 
metric and nearly symmetric double barrier resonance 
supporting tunnel junctions (DBRTS). When the asym- 
metry is stronger the noise profile has a single peak form 
with peak position at Eq- dSs/dV\ea in both cases is given 
by 16a(l - a) [1/4 - a{l - a)] (in units of 4e^/?i), which 
is for symmetric coupling case (a = 1/2). 

The sum St -\- Ss determines the observed signal. Fig- 
ure Q illustrates the above discussion using the parame- 
ters £0 = 2 eV, F = 0.04 eV, and T = 10 K. Here and 
below we use a.u. (atomic units) for current, noise and 
differential noise, e/to, e^/io, and e^ /h respectively, with 
e the electron charge, ip the atomic unit of time, and 
fi = h/2'iT the Planck constant. Fig. ^ shows plots of 
the conductance and the differential noise for symmetric 
(a = 0.5) and strongly asymmetric (a = 0.1) coupling. 
In both cases dl/dV is a Lorentzian (dotted and dash- 
dotted lines respectively) , while form of dS/dV changes 
from a two-peak structure with complete noise suppres- 
sion at fiL = So for the symmetric case to a one peak 
structure in the asymmetric situation. Note that in the 
two-peak case the peak heights are different due to the 



FIG. 1: (Color online) Conductance and differential noise vs. 
applied source-drain voltage, (a) dS{u! — 0)/dV (solid and 
dashed lines) and dl /dV (dotted and dashed-dotted lines) 
for Q = 0.5 and a = 0.1 respectively, (b) Contour plot of 
dS{u = Q)/dV . See text for units. 



thermal noise contribution. In principle both the peaks 
positions in the two-peak structure and the differential 
noise at ^l — £o may be used as sources of information 
on the asymmetry (expressed by the parameter a) in the 
molecule-leads couplings. A single peak dS/dV shape 
indicates highly asymmetric coupling. Consequently, in 
STM experiments the shape of the differential noise vs. 
applied voltage plot will change when bringing the tip 
closer to the sample, decreasing the Tl/Tr asymmetry. 

The effect of temperature is demonstrated in Fig. [21 
where results for two temperatures, T = 10 K and 
T = 100 K respectively, are presented for the symmet- 
ric coupling case, a = 0.5. The differential conductance 
peak, as is expected, becomes wider and lower with in- 
creasing T (Fig. |21i). At the same time the two peak 
structure of the differential noise plot becomes broader, 
and the difference in the peak heights of the two-peak 
structure in the differential noise plot becomes more pro- 
nounced (Fig. [33). 



eV (eV) 
4.0 




-0.1 0.1 
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FIG. 2; (Color online) Conductance (a) and differential noise 
(b) vs. applied source-drain voltage. Shown are results for 
T = 10 K (solid line) and T = 100 K (dashed line). 



WEAK ELECTRON-PHONON COUPLING 



We now turn to the model, Eq. (|ll|l . that includes 
the effect of electron-phonon interactions on the bridge. 
In this section we consider weak electron-phonon cou- 
phng which can be treated with the standard diagram- 
matic theory on the Keldysh contour. The lowest non- 
zero diagrams lead to the Born approximation (BA), 
while dressing of these diagrams yields its self-consistent 
version (SCBA). Detailed description of this approach 
and the derivation of the self-energies expressions can be 
found in our previous publication^ Within the Born ap- 
proximation the lesser and greater electron self-energies 
due to coupling to the local phonon can be expressed 
for the model at hand in terms of the (molecule pro- 
jected) electron and phonon densities of states, Pei{E) = 
-lm[G''{E)]/TT and pphiuj) = -sgn(w)Ini[i:>'-(^)]/7r re- 
spectively, and of the electron and phonon energy- 
resolved occupations, n{E) and N{uj) respectively)^^ as 



follows (for details see Ref. |5 



^phiE) 



-27riMjj / duj pph{uj) 
Jo 

X {[1 + N{Lu)]p,i{E -uj)[l- n{E - Lu)] 
+N{uj)p,i{E + uj)[l~ n{E + Lu)]} 



(27) 



S<,(i?) 



2TriM. 



-f2 



dujpphiuj) 



X {N{Lu)pei{E - uj)n{E - uj) (28) 

+ [1 + N{uj)]pei{E + Lo)n{E + w)} 

Note that in the 1-level bridge model employed here all 
the Green functions and self-energies are scalars defined 
in the bridge subspace. In the absence of electron-phonon 
coupling the densities of states are 



Pel{E) - 
Pph{^) = 



r/27r 



{E-e^Y + {T/2Y 
[lo - cJo)2 + (7ph/2)^ 



(29) 
(30) 



where F = F^ + Ffl, see Eq. (I15|l . is resonant level escape 
rate due to coupling to the leads and 7^/1 is the local 
phonon decay rate due to coupling to the thermal bath. 
The phonon occupation N{uj) in this limit is given by the 
Bose-Einstein distribution, N{uj) = ^e^^"^ — l] , while 
the electron occupation in the junction is 



n{E) = ^hiE) 



r 



■fB.{E) 



(31) 



where JkIE) is Fermi-Dirac distribution of the contact 
K. The electron relaxation rate due to coupling to the 
primary phonon is given by 



TpH{E) = i[T>^{E) 



^^h{E)] 



(32) 



The retarded electron self-energy due to this coupling 
is taken in what follows to be purely imaginary (tacitly 
assuming that its real part - level shift - has been incor- 
porated into eo), given by 



Kh{E) 



Xph{E) 



(33) 



Keeping terms up to the second order in the electron- 
phonon coupling we get for the electron Green functions 



G^{E) = GliE) ^Gl{E)^-^^Gl{E) 
G<->{E) = \Gl{E)\'{[l-iTTp,,{E)p,j{E)] 



(34) 
(35) 



'{E) 



'{E)] 



-•ph 



-{E)} 



where 



Gl{E)^[E-eo + iT/2Y^ 



(36) 



and where the lesser and greater electron self-energies due 
to coupling to contacts, S^'^ {E), K = L,R are given by 
Eqs. (P I and fl ^ . 

Using l(T^ . lfn|) . and j5^ in Q the average steady- 
state current through the junction is obtained in the form 



/ = 



2e 



-\-oc 



clE 

2^ 



T{E)[fL{E)-fR{E)] (37) 



where 

T{E)^To{E) 



{E - eof - (r/2)2 TpniE) 



iE-eor + iT/2y r 



(38) 



with To{E) given in H23|l . The first term in parentheses in 
(|38|l yields usual Landauer expression, while the second 
gives the second order (in Ma) correction associated with 
the electron-phonon coupling. As discussed in Ref. ISTl 
this correction results from both the renormalization of 
the elastic transmission channel and the opening of the 
inelastic channel. 



Using lO, linil, (ESI), ESI in (CDl) we get for the zero 
frequency noise 



where 
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One sees that in the presence of electron-phonon interac- 
tion the noise expression can not be separated into ther- 
mal and shot noise parts as is given by (|20|I - H22|) . 

Eqs. (|37|l and 1)39(1 are used in the calculations pre- 
sented below. We will consider two limits of the electron 
tunneling process: 1. resonant, F 3> |i? — £o|, and 2. off- 
resonant (superexchange), F <C \E — eo\, where E is the 
energy of the tunneling electron. The first situation was 
considered in Ref. ■55'. The second is a common situa- 
tion in molecular junctions where the HOMO-LUMO gap 
is larger than both the orbital broadening and applied 
voltage. It was pointed out by us and others that since 

(gZe°)^+(r/2)'i ^^ negative (positive) in the resonant (non- 
resonant) situations the current change upon the onset 
of inelastic tunneling, 51 ^ is negative in the resonant case 
and positive in the off-resonance limit. 

Figure El presents the zero frequency noise S{uj ~ 0), 
displayed against the applied voltage for the resonant 
tunneling case in (a) symmetric coupling, a = 0.5, and 
(b) strongly asymmetric coupling, a = 0.01, situations. 
Parameters of the calculation are T = 10 K, eg = 0.05 eV, 
F == 0.5 eV, 77 = 0.5, Sf = eV, ujq = 0.1 eV, 
-fph = 0.01 eV, and A4 = 0.1 eV. The full and dashed 
lines in Fig. El correspond to the calculated noise in the 
presence and absence of electron-phonon interaction, re- 



spectively. The dotted line results from an attempted 
approximation in which the noise was computed from 
Eqs. (EOJ-lElll except that Ta{E) was replaced by T[E), 
Eq. lEHJ- Two points are noteworthy here. First, opening 
of the inelastic channel, at V/ojq = 1, may lead both to in- 
crease (symmetric coupling, Fig.E^) and decrease (asym- 
metric coupling. Fig. EI5) of the noise signal. Second, 
modifying Eqs. ((20|I - H22|I by replacing the bare trans- 
mission Tq{E) by the phonon-renormalized transmission 
T{E) shows a similar qualitative dependence of the noise 
on the electron-phonon interaction as our more rigorous 
results, but fails quantitatively (compare full and dotted 
lines in Fig. E^). 



While the calculations presented below are made with 
no additional assumptions, a valuable insight can be 
gained by considering the low temperature limit, T -^ 
0, assuming luq <^ AE and/or wq ^ F, and taking 
the voltage large enough so that the inelastic chan- 
nel is open, /i^ — fJ-R. = |eV| > hujQ. In this case 
[fUE)-fRm'^h{E)[l-fRiE)] « \h{E)~fR{E)\ 
and /if (i5)[l — fxiE ± ojq)] ~ 0. Under these additional 
assumptions we obtain for the noise change upon the on- 
set of inelastic tunneling (up to M^) from (|39|) in the 
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FIG. 3: (Color online) Zero frequency noise with (solid line) 
and without (dashed line) electron-phonon interaction for (a) 
symmetric, a = 0.5, and (b) asymmetric, a = 0.01 molecule- 
leads coupling. Also shown is approximate result (dotted 
line) where noise is calculated as a sum of thermal and shot 
noise, Eqs. 1201 - 122(1 . with renormalized transmission coeffi- 
cient T, Eq. I38II . replacing the elastic tunneling transmission 
To, Eq. (1231 . See text for parameters. 



resonant tunneling situation 



bS(uj = 0) 



4e^ 






lUiE) - fRiE)] 2iUE) - 1) + r/i + r/|] 



To{E) 



^TlTr 
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(43) 
,Il£_r 

(44) 



Using H44|l in (|43|l one can show that for symmetric cou- 
pling, a = Tl/T = 0.5, 6S > 0. This together with 
51 < in the resonant tunneling case leads to (predicted 
in Ref. 55) increase in Fano factor, F = S/ 1, upon open- 
ing of the inelastic channel. The situation can be different 
however for asymmetric junctions. Indeed, Eqs. H43|l and 
(|44|l imply that if the inequality 



a 



a - 



9 

Vl 



Vr 



>0 



(45) 



is satisfied the change in the zero frequency shot noise 



FIG. 4: (Color online) Ratios of Fano factors (solid line), zero 
frequency noises (dotted line) , and currents (dashed line) with 
and without electron-phonon coupling vs. applied voltage 
in the resonant tunneling regime. Shown are (a) symmetric 
a = 0.5 and (b) asymmetric a — 0.01 coupling case. See text 
for parameters. 



is negative, SS < 0. In this case the Fano factor can 
decrease with the opening of inelastic channel. 

Figures 0^ and 2Jd (which were obtained, as already 
noted, without invoking these simplifying assumptions) 
validate these observations. Fig. 0Jl shows, for a sym- 
metric junction a trend (smaller current, larger noise 
and consequently larger Fano factor in the presence of 
electron-phonon interaction) similar to that suggested by 
Chen and Di VentraS (see Fig. 2 there). However the 
trend shown in Fig.^ for the asymmetric junction is op- 
posite. It should also be noticed that while thephonon 
induced change in the noise was obtained in |5a to be 
of order ~ M^, Eq. (|43|l shows contributions of order 
~ M^ that dominate in the weak electron-phonon cou- 
pling limit. 

Figure |5Ji shows change in the noise properties with 
changing asymmetry in coupling to the leads as expressed 
by the parameter a. Note that the current and noise are 
slightly asymmetric about the symmetric [a = 0.5) point. 
This asymmetry can be rationalized by noting that the 
difference F^ ^ Tr implies different tunneling distances 
between the molecule and the left and right leads, re- 







FIG. 5: (Color online) (a) Ratios of Fano factors {F/Fo, solid 
line), zero frequency noises {S/So, dotted line), and currents 
(7//o, dashed line) with and without electron-phonon cou- 
pling vs. asymmetry in coupling to the leads in the resonant 
tunneling regime, (b) F/Fo a function of two asymmetry pa- 
rameters, a and tjl in the resonant tunneling regime. Dotted 
line represents F/Fo = 1 case. Parameters are the same as in 
Fig. Eland the voltage bias is 1/ = 0.15 V [eV = 1.5a;o). 



suiting in different tunneling probabilities towards these 
leads for an electron that lost energy to molecular vibra- 
tions. Figure 03 gives a broader view of this issue. It 
shows a two-dimensional plot of the ratio of Fano factor 
with and without electron-phonon coupling, F/Fq, as a 
function of two asymmetry parameters, a and tjl- At and 
about the symmetric point a = 0.5 this ratio is greater 
than 1, but it becomes smaller (and can be smaller than 
1) in strongly asymmetric situations, as shown. 

The qualitative observations of Ref.jS^are seen to hold 
only for symmetric molecule-lead coupling but not for 
highly asymmetric situations as encountered in STM con- 
figurations. 

Consider next the off-resonant limit, F <C \E — EqI- 
In this case SI is positive. When also cjq ^ l-B — eo| 
and T — » we can use the argument above Eq. l|^ to 
show that SS > 0. In the absence of electron-phonon 
interaction Iq ^ Tq eV and ^o ~ To[l — Tq] eV one gets 
Fq ~ Sq/Iq ^ 1 — Tq. When this interaction is present 
/ ~ ro[l + rj,^/F] eV and S - {To[l-To]-f ToFp;,/F[r,i-|- 
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FIG. 6: (Color online) Ratios of Fano factors (solid line), zero 
frequency noises (dotted line) , and currents (dashed line) with 
and without electron-phonon coupling vs. applied voltage in 
the off-resonant tunneling regime. See text for parameters. 



?7|-2To]}ey, so that 
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^i 



4] <i 



(46) 



Figure El (which again does not rely on the extreme limit 
used to develop our intuition) illustrates this situation. 
Parameters of the calculation are those of Fig.^Ji except 
that £0 = 5 eV. With opening of the inelastic channel 
both current and noise grow, while the Fano factor de- 
creases. Contrary to the situation of resonant tunneling, 
asymmetry in the coupling to the leads does not provide 
qualitative difference in this case. 

Figurc[7|shows the effect on the current and noise prop- 
erties of transition from the resonant to the off-resonant 
tunneling regime as may be observed by applying a suit- 
able gate potential. Parameters of the calculation are 
those of Fig. ^ with position of the molecular orbital 
changing relative to the Fermi energy. Shown are the 
ratios of the current (dashed line), noise (dotted line) 
and Fano factor (solid line) to their values in the absence 
of electron-phonon interaction. As noticed above and as 
implied by Eq. H38|l the current goes through a minimum 
at resonance. In contrast the noise (and consequently 
the Fano factor) maximize at that point. Note that the 
results are obtained using eV = 0.15 eV > Tiwo. For 
eV < hujQ phonon effects are negligible and the resulting 
lines will approach unity. 
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FIG. 7: (Color online) Ratios of Fano factors (solid line), 
zero frequency noises (dotted line) , and currents (dashed line) 
with and without electron-phonon coupling vs. position of the 
molecular orbital relative to the Fermi energy. See text for 
parameters. 



VI. STRONG ELECTRON-PHONON 
COUPLING 

In this section we sketch the procedure and present 
results of our numerical calculations for the case of in- 
termediate to strong electron-phonon coupling. Details 
of the calculation procedure can be found in Ref. |53 and 
here we mention only several important points. A small 
polaron (canonical or Lang-Firsov) transformationS^iZ^ 
applied to Hamiltonian Hll|l leads to 



st; 



fee{L,_R} 



Y, {VkclcXa+h. 
ke{L,R} 



ujoa^a + ^ ujphybp + ^ UpQaQi) 



(47) 



where 



Af2 
eo=eo-A A « ^^; 

UJQ 



(48) 



A is the electron level shift due to coupling to the primary 
phonon and 



Xa = exp 



^0 



(49) 



is the primary phonon shift generator with Pa = 
—i (a — d^) being the primary phonon momentum. We 
are looking for the single electron Green function, ap- 



proximately given (on the Keldysh contour) by 

G'(ri,T2) = -J<T,c(Ti)ct(T2)>H 

- -r < T,c(ri)X,(Ti)ct(r2)ll(r2) >h 



(50) 



« -- < rcc(ri)ct(r2) >H < TrMri)Xl{T2) >n 

= Ge(Ti,T2) < T,Xa{Ti)Xi{T2) >h 

In Eq. H50|l Tc is the time ordering operator on the 

Keldysh contour, H is the transformed Hamiltonian, 
Eq. H47|) . and < ... >^ implies that the indicated time 
evolutions are to be carried out with this Hamiltonian. 
Also 



<T,XaiTl)XliT2)>-- 



(51) 



exp |a2 [inDp^p^{Ti,T2)- < PI >] } 



is an approximate (second order cumulant expansion) ex- 
pression for the shift generator correlation function in 
terms of the primary phonon Green function 

Dp^pSti,T2) = -J < TcPa{Ti)Pa{T2) > (52) 

n 

The (approximate) electron Green function G is ob- 
tained from a self-consistent solution of the approximate 
(second order in the coupling to the contacts) equations 
for the phonon and electron Green functions 

Dp.pAr,r')^D^^}p^{T,T') (53) 

+ fdTl I dr2D^^}pjT,Tl)Ilp^p^{Tl,T2) Dp^p^{t2,t') 

G,{T,T')^G^^Hr,r') (54) 

+ E fdri fdT2G'fHT,Tl)j:,,K{n,T2)Gc{T2,T') 

K={L,R} 

These equations (for derivation ses^) are analogs of the 
usual Dyson equation. The functions Hp^p^ and Sc,i<-, 
which are analogs of usual phonon and electron self- 
energies, satisfy the equations 

np^pAruT2) = J2\Up\'Dp,P,{rur2)-iXl Y. I^'^l' 

/3 ke{L,R} 

X %fe(r2,Ti)G,(Ti,r2) < T,Xa{Ti)Xl{T2) > (55) 

+ (ti ^ T2)] 
Sc,K(ri,T2) = Y l^fclV(n,r2) < T,XaiT2)XliTi) > 

(56) 



keK 



Here K — L,R and gk is the free electron Green function 
for state k in the contacts. 

The Langreth projections^^ of Eqs. (|53|I - H5()|) on the 
real time axis are solved numerically by using grids in 
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time and energy spaces, repeatedly switching between 
these spaces as described in Ref. [5^ After convergence, 
the Green function and self-energies needed in the noise 
expression UlUfl are available as numerical functions of E. 

In the presence of electron-phonon interactions the ex- 
pression for the average steady state current can still be 
written in a form similar to H19|! ^^ with a rcnormalized 
transmission coefficient T{E) replacing Tq{E). T{E) can 
be written as T{E) = TLTR2npei{E)/r, where Pei{E) is 
now the interacting electronic density of states.-^^ Again, 
the expression for the noise can no longer be cast in the 
simple additive structure of Eq. (pn|l iSi 

In what follows we present and discuss some numer- 
ical results based on this procedure. These results will 
serve to illustrate the sensitivity of the current and noise 
calculations to the level of approximation used, the man- 
ifestation of the phonon sideband structure in the noise 
spectrum and the conditions for observing a single or split 
central peak in the noise spectrum in strong electron- 
phonon coupling situations. We focus on these issues 
while limiting the size of our parameter space by taking 
symmetric junctions with F^ = Tn (i.e. a = 0.5) and 
rji ~ —rjn = 0.5, eo — 2 eV above the unbiased Fermi 
energies, T = 10 K, ujq = 0.2 eV, Ma = 0.2 eV and 
jph = 0.01 cV. 

FigurelHlshows the differential conductance dl/dV and 
noise spectrum dS/dV evaluated at ^ = 0.01 V and plot- 
ted against the resonance energy eg (that is in principle 
controllable by an imposed gate potential). The dotted 
lines in Fig.|Slrepresent the current and noise obtained in 
the purely elastic case (no electron-phonon interaction). 
The solid lines are results of the converged self-consistent 
procedure described above while the dashed lines re- 
sult from the lowest order approximation to Eq. H50(l 
in which the two factors (electronic and phononic) that 
enter the Green function G are calculated in the ab- 
sence of electron-phonon interaction.^^ Two important 
features should be noticed. First, because of its split 
peak structure, the noise spectrum is far more sensitive to 
the electron-phonon coupling than the conduction plot. 
Second, for the small source-drain potential used here 
no phonon sidebands appear. A similar observation was 
made before for the conduction spectrumi^SiSi This ob- 
servation stands in contrast to the results of Zhu and 
Balatsky^^ who obtained phonon sideband structure in 
the noise spectrum plotted against the gate potential. 
Indeed the procedure of Ref. ^ would lead to sideband 
structure also in the dl/dV vs. gate potential spec- 
trum and results from the semiclassical approximation 
discussed abovei^ It should be noted that this approx- 
imation is valid only when hwo ^ ksT, whereupon the 
sideband structure will be suppressed by thermal broad- 
ening associated with the width of the Fermi function. 

In contrast to gate- voltage plots made at small bias po- 
tentials, the phonon sideband structure is observed when 
the differential conductance is plotted against the source- 
drain voltage Vm^^ A similar behavior is observed for 
the noise spectrum. Figure |51 presents results of the 
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FIG. 8: (Color online) Conductance (a) and differential noise 
(b) vs. level position (gate voltage experiment). Shown are 
results of self-consistent calculation (solid line) and the zero 
order approximation (dashed line). Results in the absence of 
electron-phonon coupling (dotted hue) are given for compar- 
ison. 



corresponding calculation for the case of relatively weak 
electron-phonon coupling. Parameters of the calculation 
are the same as in Fig. |S1 except F^ — F^ = 0.08 eV 
and Ma = 0.08 eV. In this weak coupling case only one 
sideband appears, indicating a single phonon creation by 
the tunneling electron. Also in this limit the low order 
calculation appears to be sufficient — self-consistent cor- 
rections are seen to be small. 

A stronger coupling case is presented in Figure ^1 
where the parameters of calculation are the same as in 
FiglHl Several phonon sidebands indicating phonon emis- 
sion are observed on the right of the central (elastic) 
peak. The satellite feature that appears to the left of 
the central peak is obtained only in the self-consistent 
calculation (solid line). As was discussed in |5^ this fea- 
ture indicates phonon absorption by the tunneling elec- 
tron and results from heating of the junction by electron 
flux. We see again that the shape of the differential noise 
curve. Fig. llOb . is more sensitive to interaction renor- 
malization (self-consistency of calculation) than that of 
the conductance. Fig llOb . especially in its central peak 
structure. 
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FIG. 9; (Color online) Conductance (a) and differential noise 
(b) vs. applied source-drain voltage for the case of relatively 
weak electron-phonon coupling. The different lines represent 
results for the full and zero order calculation and for the case 
with no electron-phonon coupling (line notation as in Fig.jHJl. 



FIG. 10: (Color online) Conductance (a) and differential noise 
(b) vs. applied source-drain voltage for the case of relatively 
strong electron-phonon coupling. Line notation is as in Figs.|Hl 
andU 



A significant difference between the split peak struc- 
ture of the elastic central peak in Figures 1^ and [TUl (solid 
lines) and the elastic result (Fig. ^ and dotted lines 
of Figs. El and I10|l when considered in the symmetric 
molecule-leads coupling case, a — 0.5, is the fact that 
the differential noise spectrum vanishes between the two 
peaks in the absence of electron-phonon coupling but re- 
mains finite when this coupling is present. To examine 
the significance of this difference we focus on the T = 
limit and the zero order approximation. The lesser and 
greater Green functions are given by 



The retarded Green function can be estimated from 
the Lehmann representation^ 



G''(E) 



+°^ dE' G>{E')-G<{E') 
_^, ~^ E- E' + iS 



(59) 



G<{E) 



le 



-A^ Y^ 



Z — ^ 7^1 



(57) 



T^lOJPl - E - ncop) + Tr0{hr - E - najp) 
{E-ep + nwo)2 + (r/2)2 



G>{E) = -^e-^^l^^ (58) 

n=Q 

TlO{E - nujo - (II) + ^rO{E - hloq - ^ir) 



{E-eo- nujor + (r/2)2 



For the case of relatively weak coupling to the leads 
r <^ cjQ (which is the case where phonon sidebands may 
be observed), the part of C responsible for the cen- 
tral peak region can be approximated by using in H59|) 
only the n — terms of ()57|l and (|58|l . This leads to 
G';^^q{E) = e-^»[£;-eo + «r/2]-i. Using this expression 
for G^ together with G^ and G^ approximated again by 
their n = terms in Eq. 1)1 0|l leads to 
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dSn^oiuj = 0) 
dV 



4e3 



{l-6)x To(A*fl)e-^» (tjI + rjl - 27jLme'^' ~ T^km)^'^ 
I 



(60) 



where we have also used Eqs. (I13ll - (|15ll and where To is 
given by Eq. it^ . 

Eq. (inOJ is the analog of Eq. 123 (the thermal 
noise contribution vanishes in this T = limit), 
and indeed leads to H25() for \a = 0. For fi- 
nite electron-phonon coupling the result is seen to de- 
pend on the capacitance ratios tjl and rjji. The re- 
quirement that To{E)[l — To{E)] has three extrema 
for a double-peak structure to be seen in the dS/dV 
vs. V plot can now be applied for the function 

ME)e-^" {r,l +r]l- 277^77^6"^' - To(£;)e-^') . This 

leads to an analog of Eq. H26|l 



a^ - a + ^ < 



p = e^» 



(61) 

2iiLm (e^» - 1) (62) 

which reduces to (|26|l when Aa ^ 0. The peak positions 
are given by Eq T [r/2]-\/8a(l — a)/p — 1. One can show 
also that 



dS_ 
dV 

dS_ 
dV 



peak 



271% 4 



(63) 



4e 2 

— — (5e-2^» 16a(l-a)[p/4-a(l-a)] (64) 

27rn, 



In the absence of electron-phonon coupling p — 1 and the 
expressions above take the forms shown in section Hvl for 
elastic tunneling case, where for a — 0.5 4^1 =0. This 

° ' dV Ieq 

vanishing does not occur even for a = 0.5ifp^l. Con- 
sequently, the peak positions and noise values spectrum 
near and between the peaks contains information about 
the asymmetry of the coupling to the leads, capacitance 
parameters and strength of electron-phonon interaction. 
For very strong electron-phonon interaction, Aa — > cxd, 
condition H61|l can not be satisfied for any choice of the 
junction parameters. In this case only one peak in dS/dV 
in the central region will be observed. Figure ^2 illus- 
trates this effect. The parameters used in calculation 
are the same as in Fig. ^| The calculation is done for 
symmetric coupling, a = 0.5 and rj = 0.5. This behav- 
ior is seen already in the zero order approximation as is 
shown in Fig. ^2 Fig. Illb shows evolution of the two 
peak structure into the one central peak with increase 
of electron-phonon coupling, calculated in this level of 
approximation. Fig. Illb shows result of self-consistent 
calculation. A two peak structure is observed for weak 
{Ma = 0.04 eV, solid line) and a single peak for strong 
{Ma = 0.3 eV, dashed line) electron-phonon coupling. 
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FIG. 11; (Color online) Diflterential noise vs. applied source- 
drain voltage for a = 0.5. (a) Surface plot of dS{uj — 0)/dV as 
function of source-drain voltage and electron-phonon coupling 
strength. Zero-order result, (b) dS{uj — 0)/dV for weak 
(M = 0.04 eV, solid line) and strong {M = 0.3 eV, dashed 
line) electron-phonon coupling. Self-consistent calculation. 



This spectral change from two to one peak structure 
with increasing electron-phonon coupling results from 
transient localization and increasing dephasing of the 
electron on the bridge and can be used for estimating 
the junction parameters. For example, in an STM ex- 
periment, increasing tip-substrate distance will lead (for 
not too strong electron-phonon coupling) to a change of 
the differential noise/voltage profile from a two to a one- 
peak structure as was seen in Fig. ^ (the general shape of 
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Fig. n will persist also for a non-zero Ma). Eq. H61|) writ- 
ten as an equality is a relationship between a, 77, ujq and 
Ma at this transition point. Information on ujq is easily 
obtained from the inelastic tunneling structure of either 
the average current or the noise. To get information on 
a one can use Eqs. (jSJ, IJ57JI and (|5H|l to get for T ^ 



dl 



dV 2TTh r ^ ^ n\ 

n=0 



(eV — nujo) 
+9{-eV - ULJo) 
where 

T±n{E) = 
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It is easy to see from (|HK| and (|Hn|l that 



dvl 



di 
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tiL=£o+nuJo 
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dl 
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(66) 
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In addition Eq. (|18|l provides a relationship between the 
parameters a, 6 and 77. An independent determination of 
the voltage division factor 5 will therefore determine also 
the other parameters (within our simple junction model). 
We note that information on the voltage distribution in 
a junction can be obtained directly^^ or indirectly from 
I/V plots at positive and negative bias. 



VII. CONCLUSION 

We have studied inelastic effects on the conduction and 
noise properties of molecular junctions within a simple 
model that comprises a one level bridge (arbitrarily cho- 
sen above as the molecular LUMO) between two metal 
contacts. The level is coupled to a vibrational degree of 
freedom that represents a molecular vibration, which in 
turn is coupled to a thermal phonon bath that represents 
the environment. 

The case of weak electron-phonon coupling is inves- 
tigated using standard diagrammatic technique on the 
Keldysh contour as described in |53- We find that the 
predictions of Ref. ISSl increase in zero frequency noise 
and Fano factor with opening of the inelastic channel, 
hold only in a particular range of parameters — res- 
onant tunneling with symmetric coupling to the leads. 
Different modes of behavior are found upon changes that 
can be affected by moving an STM tip or changing the 



gate potential. The crossover point between the different 
modes of behavior provides information on the strength 
and asymmetry of the molecule-leads coupling and the 
position of the molecular level relative to the Fermi en- 
ergy. 

For the moderate to strong electron-phonon coupling 
case we have implemented a recently developed self- 
consistent scheme on the Keldysh contour. ^^ We have 
studied the zero-frequency differential noise spectrum for 
resonant inelastic tunneling and have investigated the 
shape of this spectrum and its dependence on the pa- 
rameters of the junction (asymmetry in coupling to the 
leads and electron-phonon interaction strength) . We find 
that 

• As in the I/V spectrum, the differential noise spec- 
trum in resonance inelastic tunneling shows a typi- 
cal structure characterized by a central feature and 
phonon sidebands. 

• Contrary to Refi^ we find that displaying the dif- 
ferential noise dS{iu = 0)/dV against the gate po- 
tential does not allow observation of phonon side- 
bands. Rather the sideband structure can be ob- 
served in the usual source-drain measurements. 

• The central feature exhibits a crossover from dou- 
ble to one central peak structure, both for increas- 
ing asymmetry in molecule-lead coupling and for 
increasing vibronic coupling strength. This change 
of spectral shape may be used to gain information 
on the parameters that determine these properties. 
In particular, asymmetry can be controlled by the 
tip-molecule distance in STM experiments. 

It should be emphasized that as in other studies of cur- 
rent noise in nanojunctions we have focused on the elec- 
tronic current in static junction structures, here aug- 
mented by electron-phonon interactions. Other sources 
of noise, e.g. structural fluctuations that affect tunneling 
matrix elements and molecule- lead coupling may exist, 
and future studies should assess their potential contribu- 
tions. 
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